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ABSTRACT 

A number of dwarf spheroidal (dSph) galaxies are known to contain a more extended, 
metal-poor population with a flatfish velocity dispersion profile, and a more concen- 
trated, metal-rich population with a velocity dispersion declining with radius. The 
two populations can be modelled with Michie-King distribution functions (DFs) in 
the isothermal and in the sharply truncated limits, respectively. We argue that the 
truncation of the metal-rich population can be traced back to the spatial distribution 
of the star forming gas. Suppose 5 is the exponent of the first non-constant term in the 
Taylor expansion of the total potential at the center [5—1 for Navarro-Frenk- White 
or NFW halos, 5 = 2 for cored halos). Then, we show that the ratio of the half-light 

radii of the populations RyI"^/ must be smaller than the ratio of the line-of-sight 

velocity dispersions (7ios,2(-Rh,2)/o-ios,i(-Rh,i)- 

Specialising to the case of the Sculptor dSph, we develop a technique to fit simul- 
taneously both populations with Michie-King DFs. This enables us to determine the 
mass profile of the Sculptor dSph with unprecedented accuracy in the radial range 
0.2 < r < 1.2 kpc. We show that cored halo models are preferred over cusped halo 
models, with a likelihood ratio test rejecting NFW models at any significance level 
higher than 0.05%. Even more worryingly, the best-fit NFW models require concen- 
trations with c < 20, which is not in the cosmologically preferred range for dwarf 
galaxies. We conclude that the kinematics of multiple populations in dSphs provides a 
substantial new challenge for theories of galaxy formation, with the weight of available 
evidence strongly against dark matter cusps at the centre. 
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1 INTRODUCTION 

The population of dwarf spheroidal satellites (dSphs) of the 
Galaxy furnishes an unprecedented opportunity to test the 
physics underlying the most extreme dark matter dominated 
systems. In particular, the p rediction of a cuspe d density 
profile for dark matter haloes (jNavarro et al.lll996l ) has nat- 
urally focused attention on the structure of the innermost 
parts. However, even though reasonable photometric and 
kinematic profiles of single stellar populations in several of 
the most luminous dSphs are available, the only clear conclu- 
sion is that the methods of stellar dynamics have so far been 
unable to prove whether t he central regions of such haloes 
are cored or cusped (e.g.. iKlevna et al.l |2002|: iKoch et al] 
I2OO7I : IWalker et aLll2009l : lAmorisco fc Evansll201ll V 

The last decade has also seen the growing real- 



(e.g.jlMaiewski et al. 1999; Klcvna et al. 


2004: Tolstov et al. 


2004 iKoch et all I2OO&: Battaglia et al. 


.2008a: lUral et al. 


2O10f). Carina. Fornax, Sculptor and Sextans all disDlav evi- 



dence for the co-existence of at least two stellar populations 
with different metallicities - usually an older, more spatially 
extended and metal-poor population and a younger, more 
spatially concentrated, metal-rich population. In all cases 
save Carina, the different stellar populations clearly have 
different kinematics: while the older component displays a 
fairly flattish velocity dispersion profile, the dispersion pro- 
file of the younger component often shows signs of a quite 
sharp decline with radius. As a consequence, the metal-rich 
population is concentrated in the central regions of the sys- 
tem only. This picture is shared by both dSphs - Sculptor 
and Sextans - where such kinematic profiles have been com- 
pletely disentangled (|Battaglia et al.H2008al . l201ll '). 
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This is consistent with a picture in which the dSph has 
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undergone mainly two different star formation epocfis, witli 
an underlying evolution driven by tidal stirring in the po- 
tential of the Galaxy. The metal-poor population may have 
undergone thermalization with the help of impulsive tidal 
stirring near the pericenters of orbits and of dynamical in- 
stabilities driven by the tidal field near the apocenters (e.g., 
iRead et al.ll2006l : iKlimentowski et alllioOTl ). This can lead 
to a nearly isothermal distribution of velocities. By con- 
trast, the metal- rich population may originate from a cen- 
trally concentrated distribution of gas, with the gas in the 
outskirts of the dSph being removed during the interval be- 
tween the two star formation events. The spatial distribu- 
tion of the gas may enforce an approximate cut in energy in 
the phase space distribution function of the younger stellar 
population. 

The coexistence of two different stellar populations in 
the same potential provides a stronger instrument to probe 
the properties of the dark mat ter halo in dSphs. A fi rst 
approach is already provided bv lBattaglia et ahl (|2008al ) in 
their analysis of the two populations in the Sculptor dSph. 
They found that, though a Navarro-Frenk- White (NFW) 
halo is still statistically compatible with the kinematical 
properties of the two stellar populations of the Sculptor 
dSph, a cored halo is a better fit in the context of the Jeans 
equations. Note however th at the theorem on th e central ve- 
locity dispersion proved by lAn fc Evanj (|2009l ) shows that 
an isotropic cored stellar density profile cannot be embedded 
in a cusped dark matter halo, so that Battaglia et al.'s so- 
lutions to the spherical Jeans equations are not all physical 

i see a lso the recent re sult provided by giotti fc Morgantil 
2010l )l. Very recently, IWalker fc Pefiarru bia (201l|) have 
used an ingenious argument to conclude that the popula- 
tions in both the Fornax and Sculptor dSphs strongly dis- 
favour cusped dark matter haloes. However, their result is 
based on th e assumption that t he mass e stimator formula 
presented in I Walker erall (|2009f ) (see also lWolf et~al]|2010l l 
can be applied separately to both metal-poor and metal-rich 
stellar populations in a dSph, a hypothesis which remains 
to be demonstrated. Although the velocity dispersion pro- 
file of the metal-poor population is consistent with a quasi- 
isothermal system, and thus wit h such a formula, that of the 
metal-rich is probably not (c.f.. iBattaelia et al.ll2008al ). 

At the price of assuming a functional form for the dis- 
tribution function of the stellar populations, we can signif- 
icantly strengthen the analysis by requiring the models to 
reproduce at the same time the luminosity and kinematic 
profiles for both populations. In doing this we aim to mea- 
sure a mass profile (at least in some radial interval), rather 
than a single t otal mass measure at some radius. B oth Jeans 
analyses (e.g. JWalker et al.ll2009l: IWolf et al.ll2010l') and sim- 
ple phase space modelling ( Amorisco fc Evansll201ll ') are able 
to measure the total dark halo mass inside some specific ra- 
dius, typically related to the half-light radius of the stellar 
population, but have failed to reconstruct the entire mass 
profile. This is because, though they are able to link the two 
dimensional scales of the dark matter density profile (the 
scalelength ro and the characteristic density po for exam- 
ple), they are unable to fix either of these numbers. This 
is very typical when the velocity dispersion profile only is 
fitted, as it is naturally degenerate in these two dimensional 
scales. However, the surface brightness profile can put inde- 
pendent constraints on the scalelength of the halo. 



We re-analyze the data for the two stellar populations 
in the Sculptor dSph by using phase space models, instead of 
the less complete Jeans analysis, adopting the general fam- 
ily of Michie-King distribution functions (King 1962, Michie 
1963). This represents a fairly flexible family of centrally 
isotropic distribution functions with an adjustable radial 
anisotropy in the outer regions. The two stellar populations 
are embedded in a dark halo whose density proflle can be 
shaped at will, cored or cusped. 

In Section [51 we introduce the observational data; in 
Section [3] we describe the motivations and the basic details 
of the models we develop in the paper. Section [3] sets out 
the properties of our models and establishes a compatibility 
criterion between the observed properties of a two popula- 
tions dSph and the density profile of the dark matter halo. 
Finally, the results are presented in Section [5] 



2 SCULPTOR DATA 

Although a number of our results on multiple populations 
in dSphs are general, it is nonetheless useful to have a par- 
ticular example in mind. For this, we choose the Sculptor 
dSph which is located at a distance D = 79.3 kpc from the 
Sun. iMaiewski et~aLl (| 19991 ) already provided evidence for 
a bimodality in Sculptor's metallicity distribution based on 
a properties of the colour-magnitude diagram. This result 
has been revis ed and extended by a number of later inves- 
tigations (e.g., Babusiaux et al.l 20051 : iBattaglia et al.|[2006l . 
l2008a l: 'Walker fc Pefiarrubia I bOllf T 

IBattaglia ot al. (2008b) used the FLAMES spectro- 
scope on the VLT in the Ca II triplet region to obtain 
high quality radial velocities for Sculptor members. To di- 
vide the stars into met al-rich and metal-poor populations, 
iBattagUa et~all l|2008al ) assumed a hard cut in metallicity. 
The metal-rich population has [Fe/H] > —1.5 and traces the 
red horizontal branch stars. The metal-poor population has 
[Fe/H] < —1.7 and traces the blue horizontal branch stars. 
Those stars in the intervening band of metallicities are not 
classified. The line of sight velocity dispersion profiles as- 
sociated with the red and blue horizontal branch popula- 
tions are displayed in Fig. [T] Battaglia at al.'s procedure 
is open to questioning, as the metallicity is correlated with 
the kinematics, and the procedure may introduce biases that 
are not easy to quantify. We explored this issue by reanal- 
izing the entire radial velocity detaset, looking for corre- 
lations between the position of the hard cut in metallicity 
and the derived kinematics. We find that the kinematics of 
the metal-poor population is robust against changes in the 
cut position, as well as against foreground contamination at 
large radii (different ways to estimate the numb er of con- 
taminants - see eqn (1) in IBattaglia et al] l|2008bl ) - do not 
result in significant variations of the velocity disperion pro- 
file). The metal-rich population is more prone to variations, 
though we find that binning size is more of an issue when 
compared to the position of the metallicity cut itself, which 
leaves the main features of the velocity dispersion profile 
unchanged if varied within 0.1 in [Fe/H]. 

Battaglia et al.l (l2008al ) also provide surface brightness 
profiles for the red and blue horizontal branch stars, which 
are also shown in Fig. [1] A photometric profile is usu- 
ally characterised by the half-light radius, namely the (pro- 
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Figure 1. Surface brightness profiles (stars arcmin ^) and projected velocity dispersion profiles for the metal-poor component (in blue), 
metal-rich component (in red) and the red giant branch stars (in black). Data from Battaglia et al. (2008a). 



jected) radius the encloses half the light. Battaglia et al. 
find the metal-poor population is well fit by a Plummer pro- 
file with half-light radius 15.1 ± 0.5 arcmin (348 ± 12 pc), 
whilst the metal-rich prefers a slightly cuspier Sersic pro- 
file with half-light radius 8.6 ± 1 arcmin (198 ± 25 pc). We 
re-compute the two half light radii directly from the sur- 
face brightness profiles in Fig. [1] using a model-independent 
technique, which does not require the assumption of any 
precise functional form for the surface brightness profiles 
themselves. We assume a Gaussian distribution for the ob- 
servational uncertainties of the surface brightness profiles 
and use them to generate syntetic surface brightness pro- 
files through a Monte Carlo technique. By direct integration 
of these profiles we infer the half-light radii of the two stel- 
lar populations as well as their uncertainties (variations due 
to the order of the interpolation of the syntetic profiles are 
accounted for in the final uncertainty itself). We obtain: 

i?h,i = 341 ± 9 pc, i?h,2 = 217± 15 pc, (1) 

for the half-light radii of the metal-poor and metal-rich pop- 
ulations respectively. 



As a flexible family of dark matter haloes, we use 

pM = \ (^) (^+(^) ) (2) 

Here, ro and po are the dimensional scalelength and charac- 
teristic density of the halo. In our models, the tidal radius rt 
of the dark matter halo is the same as that of the metal-poor 
stellar population. Even if, in nature, the truncation radius 
of the dark matter halo is larger than the stellar tidal radius, 
the potential experienced by the stars remains unchanged. 
The metal-rich population is allowed to have a smaller tidal 
radius. 

Throughout the paper, we will concentrate on the fol- 
lowing two choices: 

(a,b,c) = (1,1,2) 
(a,fe,c) = (0,2,3) 

The first choice yields the cosmologically-motivated 
Navarro-Frenk- White (NFW) model. The second is a stan- 
dard example of cored model, with the same asymptotic 
density fall-oflf as the NFW profile (p ~ r"^). 



3 MODELS 

3.1 The Dark Halo 

Both stellar populations are assumed to be embedded in a 
massive, extended, spherical dark matter halo, the density 
profiles of which can be shaped at will. The dark halo is 
truncated at the tidal radius of the metal-poor stellar pop- 
ulation. Given the high mass-to-light ratios of dSphs, we 
assume that the entire gravitational field is determined by 
the dark matter only and that the stellar components are 
just tracers. 



3.2 The Stellar Populations 

We are inspired by the fact that, as a result of relaxation 
processes, the final distribution function of the stars in phase 
space tends to near-isothermality, at least in the very cen- 
tral parts. This idea is supported by the universal flatness 
of the velocity dispersion profiles of the dSphs at least close 
to the centres (Evans, An, Walker 2009), as well as numeri- 
cal simulations of tidal stirring of dSphs (Mayer et al. 2001). 
Fig. 1 shows that for the Sculptor dSph the metal-poor pop- 
ulation has a fiat dispersion profile out to at least ~ 500 pc. 
The metal-rich profile is perhaps more debatable, though 
consistent with a flattish profile within ~ 200 pc. 

Deviations from fiatness of the velocity dispersion pro- 
files can be ascribed to either anisotropy or tidal effects. 
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In principle, since the two populations reside in the same 
potential and experience the same tidal perturbation, they 
should share the same cut in energy too if this is tidal in ori- 
gin. However, we are motivated by the idea that the metal- 
rich population has indeed been formed by gas which was 
concentrated in the central and more protected regions of 
the dSph. Thus, its cut-off in radius may be much smaller 
than the formal tidal radius. In what follows, we therefore 
assume that the more extended metal-poor population has 
the same tidal radius as the dark halo, but allow the ra- 
dial truncation of the metal-rich population to be set by the 
data. 

The Michie-King distribution function or DF (King 
1962, Michie 1963) is a function of two integrals of the mo- 
tion, the energy E (per unit mass) and the angular momen- 
tum L (per unit mass), namely 



fMKiE,L) 



P*.o 



■ exp 



( — 



ME) , (4) 



where 



fK{E) = exp 



<i>{rt)-E 



1. 



(5) 



The gravitational potential $ is, in general, the total grav- 
itational potential of the system, but here we assume that 
the stars can be treated as simple tracers. 

The dimensional free parameters of the Michie-King DF 
are as follows: 

(1) the velocity dispersion parameter a, which, for a fixed 
gravitational potential, tunes the scalelength of the associ- 
ated luminous component; 

(2) the tidal radius n, which is the external edge of the 
system, truncated by the tidal effects of the Galaxy, or by 
the original energy distribution of the star forming gas; 

(3) the anisotropy radius ra, which is approximately the 
characteristic length scale associated with variations in the 
anisotropy profile; 

(4) the characteristic stellar density /5*,o. As long as the stars 
are just tracers in the potential of the dark halo, this param- 
eter is an overall normalization, and has no other effect on 
the model. 

To our set of dimensional parameters, we can associate 
a useful set of dimensionless parameters. We use the scale- 
length of the halo ro and its characteristic density po as 
dimensional references, and denote dimensionless quantities 
by a super-script hat, so that the dimensionless half-light 
radius, anisotropy radius and tidal radius are i?h ~ Rh/ro, 
f'a ~ r^^/ro and ft — n/ro respectively. The dimensionless 
velocity dispersion is 

where, r espectively for an NFW a nd a cored halo, we have 
(see also I Amorisco fc Evand (|201ll )) 



(6) 



2nGporo, 



2TTGporo/3. 



(7) 



A Michie-King model for a single stellar population in a 
dSph is defined by the three numbers (o", f t , fa) and two scal- 
ings set by the dark halo (ro,po)- A two-population model 
for a dSph therefore has 3-1-3-1-2 = 8 free parameters. 




Figure 2. Anisotropy profiles of Michie-King models embedded 
in an NFW halo, all having .R^ = 0.2 and = 1/2, but with dif- 
ferent truncations. The dashed profile shows the Osipkov-Merrit 
limit ITOt . 



Michie-King DFs are never tangentially anisotropic, and 
range from perfectly isotropic to highly radial. This does not 
seem unduly restrictive on examining the velocity dispersion 
profiles in Fig. 1, both of which decline with radius. This fact 
does not suggest any great preponderance of tangential bias 
in the orbits. The anisotropy profiles l3{r) = 1 — a^/{2a'^), 
(where at = ag + is the tangential velocity dispersion) 
of Michie-King DFs are always isotropic at the centre, but 
can be adjusted at larger radii through the anisotropy radius 
fa. However, Michie-King models with the same anisotropy 
radius may have completely different degrees of anisotropy 
(depending on the values of the other two free parame- 
ters). Thus, in the following, we prefer to parametrize the 
anisotropy of our models using the value of the anisotropy 
profile at the half-light radius 



(8) 



Also, in order to avoid systems subject to radial orbit insta- 
bility (e.g., Barnes, Hut & Goodman 1986), we only consider 
,9 ^ 1/2. 



3.3 The Isothermal and the Strongly Truncated 
Limits 

To understand the properties of Michie-King models, it is 
helpful to distinguish two opposite physical regimes. The 
DFs can vary from exponential in the energy, typical of 
isothermal systems, to a power-law dependence in the 're- 
duced' energy _E — $(rt), typical of strongly truncated sys- 
tems. Notice that we are exploring a much wider range of 
Michie-King models than is usual in their most common ap- 
plication, the study of globular clusters. 

In the isothermal limit, most of the stars are unaffected 
by the tidal cut in King's DF ((5]), which is thus dominated 
by the exponential term: 



Mn) - ^(r)]/^^ > 1 



(9) 



An equivalent way to state this inequality is i?h ^ ft. When 
this holds true, the dependence of the Michie-King DFs on 
energy E and angular momentum L becomes close to the lin- 
ear combination E + L^ / {2r'i), which describes the Osipkov- 
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Merritt case (Osipkov 1979, Merritt 1985). Thus, in this in- 
stance, the anisotropy asymptotes to the profile 



(10) 



J.2 _|_ J.2 

This describes systems which are isotropic in the center 
(r <C Ta) and become more radially anisotropic at larger 
radii (r ^ ra). The isotropic, and therefore fully isothermal 
limit, is equivalent to additionally imposing 



-Rh/ra < 1. 



(11) 



This corresponds to the case in which the half-light radius 
of the stars is smaller than the other length-scales in the 
problem. In this regime, isothermality ensures a flat veloc- 
ity dispersion profile. Furthermore, if the stars are deeply 
embedded in their dark matter halo (i.e. <^ 1), then the 
stellar profile has a Gaussian density-law if embedded in a 
cored halo, and an exponential one if embedded in an NFW 
halo. Amorisco & Evans (2011) show that these systems sat- 
isfy: 

i?h/ro oc (aV<&o)'/* ~ [aio^V'^o)'^' , 



(12) 



where 5 is the exponent of the first non-constant term in 
the Taylor expansion of the total potential at the center. 
So, 5 = 1 for the NFW case, 5 = 2 for a cored dark halo. 

The opposite of the isothermal regime is the case in 
which the exponential term in King's DF ([5]) is of order 
unity. This limit describes strongly truncated systems, in 
which the dimensionless half-light radius Rh and the dimen- 
sionless tidal radius ft are of the same order. The asymptotic 
expression for the anisotropy is: 

2 <l>(rt) - ^(r) 



/3truiic('') 



9r2 



(13) 



This is true for any Michie-King model close to its tidal edge. 
So, any Michie-King model becomes isotropic on approach- 
ing its tidal radius. This is evident from the anisotropy pro- 
files plotted in Fig. [51 

Strongly truncated isotropic models are essentially 
scale-free models, whose properties depend only on the grav- 
itational potential <l?(r) and the tidal radius rt. The radial 
velocity dispersion is 

a^r) = ^(Hn) - ^r)) , (14) 

which generates cut-ofTs in the line of sight velocity disper- 
sion profiles. For the stellar density distribution, we have 



(15) 



For fixed half-light radius _Rh and degree of radial 
anisotropy, the tidal radius realized by strongly truncated 
systems is the smallest in the entire family of Michie-King 
models with a given gravitational potential. Thus, we use 
this limit to define the function 



rr'"{Rh,P) = Urn r\(a,i?h,^) 



(16) 



It is often convenient to express the tidal radius as multiples 
of thi minimum value defined by eqn. p6p . that is 



rt 



^ 1 



(17) 



where nt is the tidal number. This is useful in presenting 
our results, as 1/nt lies between and 1. 





Figure 3. The line of sight velocity dispersion profiles generated 
by Michie-King DFs embedded in a NFW dark matter profile 
(upper panels) and in a cored dark matter halo (lower panels). 
From left to right, ijjj = 0.1, iJh = 1. In each panel, profiles 
of different colours have different tidal radii, beginning with the 
smallest possible r™'"(iliii /9), in blue. The two profiles with the 
same colour in each panel have different radial anisotropics: ^ = 
and i3 = 0.5, at the extremes of our domain. 





Ri, = 0.1 


= 0.3 


Rh = l 


'5"los(-Rh) 


1.6% 


1.5% 


1.7% 




4.0% 


6.1% 


4.5% 


0"los(-Rh) 


2.0% 


1.5% 


1.2% 




5.9% 


4.6% 


3.7% 



Table 1. The relative uncertainty obtained using the luminosity 
weighted average or the value at ij^ in the mass estimator. In 
each instance, the upper line refers to NFW haloes, the lower to 
cored haloes. 



The Michie-King model for a single population embed- 
ded in a dark halo can be described by the natural set of 
3 parameters, namely dimensionless velocity dispersion a, 
tidal radius ft and anisotropy radius fa. An alternative, but 
equivalent set, is the dimensionless half-light radius Rh, the 
tidal number nt and the anisotropy parameter at the half- 
light radius /?. It is this set that we will use in the rest of 
the paper. 



4 PROPERTIES OF THE MODELS 
4.1 Dark Halo Mass Estimators 

A main aim in studying stellar populations in dSphs is to 
map the underlying dark halo. A number of authors have 
recently suggested that the mass of the dark halo within a 
characteristic radius related to the half-light radius is well- 
constrained (e.g.. Walker et al. 2009, Wolf et al. 2010, Amor- 
isco & Evans 2011). These mass estimators all take the gen- 
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Figure 4. The relation (with scatter) between the dimension- 
less half light radius R^^ and dimensionless line of sight velocity 
dispersion at R^^, <3"ios{^h) fof Michie-King models embedded in 
an NFW halo (upper panel) and in cored halo (lower panel). In 
both panels, the dashed line represents the asymptotic analytic 
relation provided by the isothermal limit. 

eral form (c.f. Illingworth 1976) 

M(AJ?h) = §i?h(aL) (18) 

where K and A are two dimensionless constants and (crj^Qj,) 
is an appropriate 'average' of the square of the line of sight 
velocity dispersion. A successful mass estimator needs to fix 
each of these three ingredients. 

Let us begin with the average (o-fos)- In the isothermal 
limit, the profile is almost flat by definition, so we do not 
need to worry at which radius we pick the line of sight veloc- 
ity dispersion. When the dispersion profile falls with radius, 
it is a moot point as to how best to provide the 'average'. 
This topic assumes some importance as it is the heart of the 
argument provided by Walker & Penarrubia (2011) in their 
analysis. 

To bring the problem of the 'average' into sharp fo- 
cus, consider two stellar systems embedded in the same 
dark matter halo and with the same half-light radius, but 
characterized by a different DF in phase space, and thus 
having different kinematics and density profiles. For exam- 
ple, consider a radially-biased and a tangentially-biased sys- 
tem. A proper recipe for an efficient mass estimator should 
pick the same 'average' velocity dispersion for both systems, 
though they will in general differ both in their velocity dis- 
persion profile aios{R) and their surface brightness profile 
E(7?). Wolf et al. (2010) suggested that a good choice is the 
luminosity-weighted average line of sight velocity dispersion 
{o'iob(R)^)t,{r)- However, a radially-biased model will have 



both higher velocity dispersion profile and higher surface 
brightness profile in the central regions and lower in the 
outer regions, as compared to the tangential model, which 
militates against the usefulness of a luminosity-weighted av- 
erage. 

Fig. [3] illustrates the dimensionless line of sight velocity 
dispersion profiles for models that span the entire physical 
range of Michie-King models, for NFW and cored halos (up- 
per and lower panels respectively). The panels show models 
with different half-light radii, tidal radii and anisotropy. In 
each panel, all the models have the same dimensionless half- 
light radius Rh and should then provide the same mass esti- 
mate and 'average'. A quick look makes it evident that the 
easiest choice is the line of sight velocity dispersion at one 
-Rh, where all the profiles roughly cross. The effectiveness of 
this choice is readily quantified. For three different values of 
the dimensionless half-light radius -Rh = 0.1, -Rh = 0.3 and 
-Rh ~ 1, Table [1] collates the relative uncertainty obtained 
by comparing the line of sight velocity dispersion at one 
half light radius (Tios(Ah) of all possible Michie-King mod- 
els (0 < /? < 0.5, < 1/nt < 1). The same uncertainty 
is calculated for the luminosity averaged velocity disper- 
sion \/ {o"ios(-R)^)s(.R) • Our recipe always performs better, 
with the luminosity- weighted average {o"1os(-R)^)e(-R) usu- 
ally showing discrepancies higher than 4%. 

Fig. |4] displays the functional relations obtained by as- 
sociating each dimensionless half-light radius -Rh to the di- 
mensionless line of sight velocity dispersion aios{Rh) for an 
NFW halo (upper panel) and for a cored halo (lower panel). 
In both panels the dashed line represents eqn. (|12p . which 
is valid in the isothermal limit for -Rh <C 1. The displayed 
shaded areas correspond to the error bars as in Table 1, cal- 
culated by comparing the line of sight velocity dispersions 
i3"ios(^?h) for models with different /? and nt, but identical 
half-light radius -Rh- The se relations arc similar to the ones 
displayed in Fig. 2 of [Amorisco & Evans (201l|), which were 
derived for isotropic King models in the isothermal regime 
only. 

The functional relations (Tios(-Rh) allow us to calculate 
directly the two dimensionless quantities A and K which 
define the mass estimator (llSp . We find 

A = 1.67 ±0.04, .fC = 5.85 ± 0.2. (19) 

This result is in agreement with I Amorisco fc Evans! l|201ll ). 
but it is worth noting its much wider range of validity, which 
covers the entire physical range of Michie-King models. Sum- 
marizing, we obtain 

M[(1.67±0.04)i?h] = (5.85^0.2) '^''""'°°^^'^^" , (20) 

which can be considered trustworthy even when radial 
anisotropy and substantial deviations from isothermality 
may be present. 



4.2 A Consistency Criterion for Two Populations 

As the two stellar populations reside in the same dark matter 
halo, their dimensionless half light radius -Rh and line of sight 
velocity dispersion o'ioa(^?h) must be picked from the same 
functional relation in Fig. [4] This means that the ratios of 
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Figure 5. In black, with uncertainty range in gray; the ra- 
tio fco- = i5'iog(iJii 2)/'''ios(^h,i) against -Rh,i when the ratio 
Rjj 1 is assumed from the available data. Upper panel: NFW 

halo; lower panel: cored halo. In both panels the yellow shaded 
area represents the poor evidence for ka provided by the observed 
line of sight velocity dispersion profiles; the dashed line represents 
the asymptotic value provided by the isothermal limit. 



(21) 



the half-light radii and the velocity dispersion 

7 _ Rl\,2 , _ '7-ioB.2(-Rh,2) 

f^R = ^ 1 i^tT = 7^ r, 

Uh_,i crios^l^Kh.lJ 

are not independent quantities. They are related to one an- 
other through the density profile of the dark matter halo 
itself. In general, we must have 

2 



Rh.2 

Rh.i 



0"los,2(.Rh,2) 



(22) 



with equality attained in the isothermal limit. This can be 
understood from Fig. 4, as the ratio fih,i/o"ios(Ah,i)^ is an in- 
creasing function of i?h,i and attains unity in the isothermal 
limit for strongly embedded systems. 

Let us now specialise to the case of Sculptor. For 
Battaglia's two-population dissection, using eqn (1), we find 
that 



kR = 0.64 ± 0.05 



(23) 



This measure (with its uncertainty) and the relations in 
Fig. |4] (with their uncertainties) generate the black line (and 
gray shaded areas) in Fig. O All the pairs {Rh,i, ka) in the 
gray areas are compatible with an NFW halo (upper panel) 
or with a cored halo (lower panel) for some value of the halo's 
scalelength ro. Also displayed in the figure as the asymptotic 
dashed line is the isothermal limit. 

Measurement of the ratio fco- remains a challenge despite 



substantial improvements in the quality of the photometric 
and spectroscopic data in recent years. We interpolate both 
the velocity dispersion profiles, and the profiles defined by 
the uncertainties, and obtain 



k^ =0.55 ±0.12 



(24) 



which is displayed in Fig. [5] as the yellow shaded area. 

It is evident that the cored halo provides a much bet- 
ter interpretation of the available data, as there is a formal 
agreement between models and data which is not the case 
for the NFW halo. Furthermore, we can extrapolate this re- 
sult and state that any other cusped dark matter density 
profile with a cusp steeper than the NFW prescription is 
going to deliver even worse results. 

In passing, we note that our consistency crite- 
rion Il22p can be reformulat ed to provide the one used by 
IWalkCT fc Penarrubia I (|201l'). In fact, if two populations fail 
to satisfy the inequality (I22p for a given S, then the loga- 
rithmic 'slope' r, estimated using eqn (llSp . will satisfy the 
inequality 



F = 



In [M(AJ?h.2)/M(Aflh.i)] 



1 + 2 



In ka 



>S+1, (25) 



In {Rh,2 / Rh.i) ' Infcfl 

thus confirming that such a system is inconsistent with the 
density profile p ~ r^~^. 



5 PHASE SPACE ANALYSIS: SCULPTOR 

To perform a complete maximum likelihood analysis of the 
Sculptor data using Michie-King DFs necessitates an explo- 
ration of an eight-dimensional parameter space. In fact, the 
main body of this analysis can be successfully reduced to the 
parallel and independent studies of the two metal-poor and 
metal-rich stellar populations, ignoring at first that they re- 
side in the same dark matter halo. The combination of the 
evidence coming from the two stellar populations then con- 
strains the range of feasible dark matter haloes. 

The details of the fitting of the Michie-King models to 
the kinematics, the photometry and - in the case of the 
metal-poor population - the tidal radius are given in Ap- 
pendix A. Here, we note that the fit to the metal-poor pop- 
ulation yields a Xs , Xa- and Xrt 1 whilst the metal-rich yields 
has only a and Xa-- This is because the tidal radius of 
the metal-rich population is free to vary, whilst that for the 
more extended metal-poor must be the same as the under- 
lying dark halo. 

The combination of all the available evidence for the 
metal-poor stellar population is provided by the likelihood 
and associated x^-squared values: 



LMpiRh,nt, l3,ro, po) — 



Xmp 



2 , 2 



Xr 



(26) 



However, for the sake of clarity, we will also consider the 
reduced likelihood 



LMp{Rh,nt,P,ro,po) = Lt. ■ La 



(27) 



partly to check that the fitting of the tidal radius does not 
have an undue and distorting effect on our results. Similarly, 
for the metal-rich population, we can define 



LMB.{Rii,nt, P,ro, po) 
Xmr 



Lt. ■ Lcr , 

2 , 2 
Xe + X<T 



(28) 
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Figure 6. The probability distributions associated with the like- 
lihoods I/MP and Lmr for the free parameters of the model when 
embedded in a NFW (left) and cored (right) halo. The blue (red) 
profiles are associated with the metal-poor (metal-rich) stellar 
component. Dashed and full blue profiles show the results deriving 
from the total likelihood II26I I and the 'reduced' likelihood II27II . 



5.1 Independent Analysis of the Two Populations 

Table [2] collects the results pertaining the metal-poor and 
metal-rich stellar populations when fitted separately, thus 
allowing each of them to have different dimensional scales 
ro and po- For each of (x|, xl + xl, xl + xl + Xrt), the 
values associated respectively with the best fit model, the 
68% and the 95% confidence regions are displayed. 

A comparison between the values for an NFW and cored 
halo makes clear that the main difference is in the reproduc- 
tion of the stellar surface brightness distribution, with the 
cored halo performing substantially better. Also, the fact 
that the increase in is very similar for an NFW and cored 
halo shows that the shape of the velocity dispersion profile 
is largely insensitive to the dark halo properties. 

For the metal-poor (metal-rich) population. Fig. [S] dis- 
plays in blue (red) the probability distributions of the free 
parameters of the Michie-King models, as deduced from 
the projection of the total likelihood. First, while the cored 
halo selects a Michie-King model with a similar scale length 
(i?h ~ 1), an NFW halo needs to be several times more ex- 















AE ^ Aa 


AE ^ Act ^ Art 


NFW 


69. i, 41.5, 45.7 


48.2, 66.9, 68.3 


49.0, 67.8, 69.7 




1.97, 2.08, 2.29 


1.12, 1.56, 1.59 


1.14, 1.58, 1.62 


cor6d 


32.7, 36.3, 39.7 


39.0, 54.3, 59.6 


40.3, 55.7, 60.5 




1.63, 1.81, 1.98 


0.91, 1.26, 1.39 


0.94, 1.29, 1.41 


NFW 


2.3, 5.6, 9.1 


3.4, 11.1, 13.8 






0.38, 0.93, 1.52 


0.34, 1.11, 1.38 




cored 


1.1, 3.3, 4.9 


2.1, 6.8, 9.9 






0.18, 0.55, 0.82 


0.21, 0.68, 0.99 





Table 2. Results of the independent analysis of the metal-poor 
(upper) and metal-rich (lower) stellar component. The table gives 
the values of the x^-quantities referring, in order, to the best fit 
models, to the 68% and to the 95% confidence regions, which have 
been calculated using the likelihoods in eqns I IA3I IA5I IA9I1 . Also 
the corresponding reduced-x^-quantities are reported in the sec- 
ond line for each halo model. Uniform priors have been assumed 
in the range defined by < Ah < 2, ^ l/nt < 1, < ^ < 0.5, 
as well as for ro and log po . 

tended {Rh < 0.2) than the distribution of stars in order to 
fit (yet with difficulties) the surface brightness profile of the 
metal-poor stars. Also, in the cored case, an almost isotropic 
metal-poor stellar population is preferred (P ~ 0.1), whilst 
an NFW halo favours a mild radial velocity bias (/? ~ 0.25). 
The metal-rich population requires a high degree of radial 
anisotropy (/3 > 0.4), irrespective of choice of halo. Note also 
the necessity of different tidal radii in the phase space distri- 
butions is supported by the results. In both cases, rit < 1.4, 
suggesting that the metal-rich stellar population is indeed 
in a different physical regime compared to the metal-poor 
one, which has nt > 3. Finally, the probability distributions 
for the dimensionless half-light radius Rh of the metal-rich 
population show that little help in determining the charac- 
teristic radius of the halo can be expected from a population 
in this heavily truncated regime. 

5.2 Two Populations in the Same Halo 

Figs [7] and [8] display the projections of our likelihoods onto 
the halo plane (ro, po), for an NFW and a cored halo respec- 
tively. In both Figs, the upper panel shows the 68 % and 95 
% confidence regions associated with the separate analysis of 
the two stellar populations with blue (red) representing the 
metal-poor (metal-rich) respectively. The full thick blue and 
red lines in the upper panels of both figures show the trends 
of the po(»'o) functional relations obtained respectively for 
the metal-poor and metal-rich stellar populations by using 
the relations iTios(i?h) displayed in Fig. |4] for an NFW and a 
cored halo. 

The lower panel shows how the independent results 
from the two stellar populations combine. We in fact im- 
pose the condition that the metal-poor and metal-rich stellar 
population live in the same dark matter halo (ro, po) via: 

itot(mi, 7712, ro, po) — iMp(Tni, ro, Po) ■ iMR(Tn2, ro, Po) 
Xtot = Xmp + Xmr , (29) 

where nii = (Rh, nt, P) is a general Michie-King model. Full 
dots indicate the best fit models. 
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R/kpc R/kpc 

Figure 9. The photometric and kinematic profiles for the two Sculptor populations in an NFW halo. Thick profiles: final fits provided 
by the best fitting Michie-King models associated with the black dot in Fig. [7] Dashed profiles: best fits provided by the Michie-King 
models, but when (r'o,po) is allowed to vary separately for the two stellar components. Surface brightness profiles are measured in stars 
arcmin"'^. 




R/kpc R/kpc 
Figure 10. As Fig.|9l but for a cored halo. 



Finally, Table 3 collects, for an NFW and cored halo, 
the values of Xtot associated with the best fitting models, 
the details of which we report here. For an NFW halo, we 
find: 

(ro,po) = (1.4 kpc,2.3 x W Mq kpc"^), 
mi = (^h = 0.256, nt = 10.0, ,9 = 0.20), 
m2 = (Ah = 0.136, nt = 1.1, = 0.37) . (30) 

In physical units, the Michie-King DF for the metal-poor 
population has a tidal radius of rt of 12.2 kpc, an anisotropy 
radius ra of 720 pc and a velocity dispersion a of 11.0 kms~^. 
The same quantities for the metal-rich population are n of 
750 pc, Ta of 140 pc and a velocity dispersion a of 15.6 
kms~^. For a cored halo, we find: 

(ro,po) = (0.36 kpc,2.1 X lO*Af0 kpc"^), 
mi = (7?h = 0.98, nt = 6.7, /? = 0.1), 
ma = (7?h = 0.54, nt = 1.1, /? = 0.48) . (31) 

This means that the Michie-King DF for the metal-poor 
population has r% — 7.8 kpc, ra — 1.0 kpc and (tq — 10.8 



kms ^. The metal-rich has n — 750 pc, = 120 pc and 
a = 16.0 kms"^. 

There are a few points that it is worth noting. As re- 
marked already, neither the velocity dispersion profile nor 
the surface brightness profile are able to constrain the char- 
acteristic scalelength of the NFW dark halo. This is rein- 
forced by Fig. [T] in which only the lower limit r^^^ > 1 kpc 
is fixed. The contours of the confidence regions associated 
with Ltot do not close as ro becomes large. The contours do 
follow accurately the functional relation Po^^ (ro) described 
by the blue or red thick full curves. By contrast, the confi- 
dence regions for the metal-poor populations do close for a 
cored halo. This is a consequence of the surface brightness 
profile, which picks out a limited interval in the scalelength 
ro, thus constraining at the same time the central density po. 
The confidence regions for the metal-rich populations still 
do not close, though this is a characteristic of the strongly 
truncated limit of the Michie-King models. 

In both Figs [7] and [S] only marginal agreement can be 
found between the 68% confidence region associated with 
the metal-poor and the metal-rich stellar populations. How- 
ever, the importance of this is minimal in the case of a cored 
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Figure 7. The 68% and 95% confidence regions in the (ro,po) 
plane for an NFW halo. The blue (red) contours are associated 
with the likelihood Lmp (^mr) for the metal-poor (metal-rich) 
stellar population. Dashed and full blue profiles show respectively 
the results deriving from the total likelihood I I26II and the reduced 
likelihood l|27|l . The black lines display the contours of the joint 
evidence. Full dots indicate the best fit models. 
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Figure 8. As Fig.[7l but for the cored halo. 



halo, since it does not cause any significant increase in the 
area of the confidence regions (or in the chi-square values) 
associated with the joint likelihood. Note, in fact, that the 
values of total Xtot associated with the joint evidence Ltot , 
do prefer a dark matter density profile with a constant den- 
sity core. 

Following lEadie et al.l (|l97ll ) , we define a set of models 
with nested parameters, by considering halos with a 1 — A 
cusp and a density fall of p{r) ~ r~'^ at large radii. This 
obviously includes the cusped NFW density profile (A = 0) 
and the cored profile (A = 1). Since the difference in best-fit 
total values between the NFW case and the cored case is 
58.5 — 46.3 — 12.2, we can reject a pure NFW density pro- 
file at any significance level higher than 0.05%. Completely 
analogous results have been obtained when the original kine- 
matic profiles taken from Battaglia et al. (2008) have been 
replaced with kinematic profiles derived directly from the ra- 
dial velocity dataset using slightly different metallicity cuts 
(within 0.1 [Fe/H]) or membership criteria (see also discus- 
sion in Sect. 2). Furthermore, although not reported here in 
detail, we investigated the issue of a milder cusp by consid- 
ering separately the A = 1/2 case. A milder cusp performs 
better than the NFW halo, but we can still reject it in favour 
of a cored halo at any significance level higher than 1.4% (we 
find a best-fit total of 52.3). All available indications sug- 
gest a monotonic behaviour of the best-fit total x^ with the 
cusp index. 

Furthermore, the best fitting NFW model (|30p . together 
with the parts of the halo (ro,po) plane supported by the 
joint evidence Ltot, sit in a region in discord with the stan- 
dard expectation of cosmological models. The concentration 
of the best fitting NFW halo is in fact as low as c ~ 17, 
whereas present-day dwarf galaxies are expected to have 
higher conce ntrations, c ~ 35 acc ording to Colin et al. 2004 
(but see also iMaccio et all (120071 ) 1. 

Figs [9] and [10] show the quality of the final fits to the 
surface brightness profiles and velocity dispersion profiles 
respectively in the case of an NFW halo and a cored halo. 
The full profiles are the best joint fits, that is obtained by 
using the same (ro,po) halo for both the metal-poor and 
the metal-rich stellar population. The dashed profiles show 
instead the best separate fits. It is evident from these plots 
that the real damage to the prospects of fitting the data 
with an NFW halo come from two sources: it is the combi- 
nation of photometry and kinematics that does the damage. 
First, although the properties of the metal-poor population 
are reproduced reasonably well, the metal-rich population 
is predicted to have a higher central surface density than 
the data. Second, the NFW halo requires a smaller differ- 
ence than the observed one between the values of the metal- 
poor and metal-rich velocity dispersion in order to repro- 
duce the observed half-light radius. This smaller difference 
is a rephrasing of the inconsistence found in Fig. 5 through 
our consistency criterion (|22p . However, note that the addi- 
tion of the information about the surface brightness profiles 
of the two stellar populations - which was not included in 
the arguments presented in Section 4 - has in fact changed 
the range of best fitting cored models. While the consistency 
criterion alone would select haloes with a larger scalelength 
{Rh ~ 0.5, see Fig. 5), the complete maximum likelihood 
analysis prefers models in which, approximately, mass fol- 
lows hght {Rh ^ 1). 
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Atot.bf 


Xtot 


NFW 


58.5 


78.1, 80.2 


cored 


46.3 


62.2, 68.0 



Table 3. The values of Xtot associated with the best fitting halo, 
as well as the 68% and 95% confidence regions, determined ac- 
cording to the likelihood 1291 1. 



6 THE MASS PROFILE FOR SCULPTOR 

Fig. [TT] shows the final mass profiles for the Sculptor dSph. 
In both upper and lower panel, the coloured shaded areas 
represent, at each radius, the 1-a uncertainty obtained for 
the mass profile using the joint evidence Ltot, for the case 
of an NFW halo (in green) and for a cored halo (in yellow). 

In the cored case, we expect to be able to measure a 
mass profile rather than a single mass enclosed within some 
radius. This is because we are able to fix the best (ro,po) 
in the halo plane. What is more surprising is that the same 
thing is true for the NFW halo, even though we have only a 
lower limit for the characteristic radius vq. However, this can 
be understood by noticing that, as anticipated in Amorisco 
& Evans (2011), the mass profile Af(r, ro) generated by the 
functional relation po('"o) becomes insensitive to the charac- 
teristic radius ro when Ri^ <^ 1. 

The points displayed in the upper panel show how our 
results compare with previous mass measures for Sculptor by 
Strigari et al. (2007, 2008), Walker et al. (2009) and Amor- 
isco & Evans (2011). The lower panel is a zoom for a smaller 
range of radii. We draw attention to the fact that, even if 
the matter of a core or cusp at the centre is considered un- 
resolved, we are still able to measure a mass profile rather 
than the single enclosed mass datapoint at a particular ra- 
dius. The mass profiles pertaining to the cored and NFW 
halo, in fact, agree within their 1-a uncertainty over the en- 
tire interval 200pc < r < 1.2kpc, thus proving a wide and 
unprecedented radial coverage. The pink sh aded area repre- 
sents t he best-fit cored halo as reported in iBattaglia et ahl 
l|2008al ). 

Also displayed in the lower panel of Fig. [TT] are 
the two mass estimates obtained using the formula sug- 
gested by Wolf et al. (2010), and applied separately to the 
metal-poor and metal-rich stellar populations. We calcu- 
late each of the two luminosity-averaged velocity dispersions 
•y/ (o"ios(-R)^)E(ii) using the observed profiles and a Monte 
Carlo technique which assumes Gaussian uncertainties. Note 
that the logarithmic slope of the mass profile we would de- 
duce from such a technique is steeper than the one we actu- 
ally measure. 



7 CONCLUSIONS 

The suggestion that stellar kinematics in dwarf spheroidals 
(dSphs) may enable the structure of the dark matter halo 
to be mapped out is at least a decade old. It motivated 
programs to gather radial velocities in the bright dSphs with 
4m and 8m class telescope and multi-object spectrographs. 
Despite a lot of work, it has proved hard to distinguish the 
structure of the dark matter halo, especially at small radii. 
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Figure 11. Upper panehThe mass profiles with errors as deduced 
from Figs. 7 and 8 for the NFW halo (in green) and the cored halo 
(in yellow). In the top panel, Strigari et al.'s (2007, 2009) mass 
measures at 100 pc, 300 pc and 600 pc are shown as black dots. 
Walker et al.'s (2009) measure at = 260 pc is shown as a blue 
dot and Amorisco & Evans (2011) measurement at 1.7i?h ~ 442 
pc is shown as a triangle. Lower panel: A zoom from the top panel 
displaying the agreement in the range 0.2 < r < 1.2 kpc, together 
with the masses estimated using Wolf et al.'s (2010) formula. 



where cold dark matter (CDM) theories predict 1/r density 
cusps (e.g., Navarro, Frenk & White 1996). 

For example, Gilmore et al. (2007) provided a succinct 
summary of the then available observational evidence. They 
noted that analyses based on the Jeans equations for six 
dSphs (Draco, Ursa Minor, Carina, Sextans, Leo I and Leo 
II) tended to suggest that the dark matter density was shal- 
lower than the 1/r density cusp. However, they acknowl- 
edged that the observed velocity dispersion profiles and 
cored fight distributions were likely to be consistent with 
dark matter halos with both central cores and cusps. 

A thorough examination of the evidence using the Jeans 
equations for the eight brightest dSphs (the six studied by 
Gilmore et al. together with Fornax and Sculptor) was car- 
ried out by Walker et al. (2009) using datasets of thousands 
of stars. They found that the most meaningful constraint 
provided by a Jeans analysis was on the total mass enclosed 
within the half-light radius (see also Wolf et al. 2010). In 
particular, the Jeans equations were not sufficiently pow- 
erful to distinguish between cored and cusped dark matter 
haloes. A companion theoretical investigation by Evans, An 
& Walker (2009) made explicit the limitations of the Jeans 
equations. They showed how assumptions as to the surface 
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brightness and profile velocity anisotropy could allow the 
misleading identification of a core or a cusp. Nonetheless, 
this is a consequence of the interlocking assumptions made 
in the Jeans analysis, and not the data, which are consistent 
with both cores and cusps. 

However, in at least two dSphs, there appears to be 
direct evidence favouring dark halos with constant den- 
sity cores. First, Kleyna et al. (2003) discovered kinemati- 
cally cold substructure in the Ursa Minor dSph. They ar- 
gued from numerical simulations that the persistence of 
such cold substructure is incompatible with the cusped pro- 
files of cold dark matter cosmogonies. They provided an 
interpretation in terms of a stellar cluster dissolving in a 
cored dark matter halo. Sec ond, Goerdt et al. (2006) and 
ISanchez-Salcedo et all (|2006i ) studied the five globular clus- 
ters present in the Fornax dSph, which orbit at projected 
radii ~ 1 kpc from the centre. In a cusped dark matter halo, 
they would sink to the centre under the action of dynami- 
cal friction in a few Gyr. In a cored dark matter halo, the 
dynamical friction timescale is much longer and so this pro- 
vides a natural resolution of the survival problem of Fornax's 
globular clusters. 

The exploration of multiple stellar populations in dwarf 
spheroidals is only just beginning. But there are now clear 
hints that it might provide a resolution of this problem. The 
fact that both populations reside in the same dark matter 
halo significantly reduces the degeneracies. To date, there 
have been three different attempts at modelling - each with 
their own strengths and limitations. It is striking that all 
three have come to the conclusion that the data are much 
more consistent with cored dark matter haloes than cusped 
one. 

Battaglia et al.'s (2008) investigation of the case of 
Sculptor used the Jeans equations and came to the con- 
clusion that the data are statistically consistent with both 
cored and cusped haloes, though the NFW profile yields a 
poorer fit for the metal-rich stars. The drawback to this 
analysis is that it is based on the Jeans equations alone, and 
there is no guarantee that the underlying mo dels exist. In 
fact, t he central velocity dispersion proved bv lAn fc Evans! 
l|2009ll shows that an isotropic cored stellar density profile 
cannot be embedded in an Navarro-Frenk- White halo, and 
so there is no solution for the phase-space distribution func- 
tion for the metal-poor population, despite the existence of 
the Jeans solution. Nonetheless, Battaglia et al.'s (2008) ba- 
sic conclusion - that a cored halo model provides a better fit 
than an NFW halo model - seems robust and is recovered 
by our analysis. 

Walker & Peiiarrubia (2011) developed a sophisti- 
cated statistical method for assigning probabilities as to 
whether dSph stars belong to metal-rich and metal-poor 
sub-populations. This is a significant improvement on the 
simple cut in metallicity used by Battaglia et al. (2008), 
which can introduce biases when metallicity is correlated 
with kinematics. As output from their statistical analysis, 
they obtain a measure of the average velocity dispersion 
of the metal-poor and metal-rich populations. This is used 
with a mass estimator of the form 1)18^ to obtain the mass 
enclosed at two points, namely the half-light radius of the 
metal-poor and the metal-rich populations. From this, the 
gradient of the dark halo mass profile is calculated. For the 
Sculptor and Fornax dSphs, they conclude that cusped NFW 



profiles are ruled out. The modelling is then simpler than a 
Jeans analysis, but Walker & Penarrubia (2011) argue that 
it is more robust. One qualm is that the method assumes 
that the two components are drawn from constant velocity 
dispersion populations, both in the statistical analysis itself 
and in the application of the mass estimator formula. Whilst 
this seems reasonable enough for the metal-poor populations 
in dSphs, it is not at all clear that the metal-rich population 
can be so described. In fact, the combined or total velocity 
dispersion of the total population must then decline out- 
wards, which is certainly not the case in, at least, Fornax. 
Nonetheless, Walker & Peiiarrubia have tested their reason- 
ing against synthetic datasets generated from a wide range 
of phase-space distribution functions, and it performs very 
well. 

Finally, in this paper, we have used the Sculptor dataset 
provided by Battaglia et al. (2008) to build phase-space dis- 
tribution functions for both the metal-poor and the metal- 
rich populations under the assumption of cored and cusped 
dark matter haloes. Although both can fit the data, cored 
dark matter haloes are preferred. There is clear evidence of 
discrepancies in the surface brightness profile of the metal- 
rich stars near the centre for cusped haloes. Notice that 
it is the combination of surface photometry with the stel- 
lar kinematics - rather than the kinematics alone - that is 
providing the decisive evidence against cusped haloes. Even 
more worryingly, the Navarro-Frenk- White models that are 
the best-fit are much less concentrated than the cosmolog- 
ical predictions from cold dark matter theories. Typically, 
concentrations of c ~ 35 are predicted from simulations 
of galaxy formation (e.g., Colin et al. 2004), as compared 
against c « 17 found in our analysis. 

The fact that all three analyses - which make different 
assumptions and have different strengths and weaknesses - 
come to the same conclusion is very striking. Further mod- 
elling of other dSphs with multiple populations is needed to 
confirm the results. However, the evidence that the dark 
haloes of the dSphs around the Milky Way do not have 
cusped form, but have cores, is now beginning to look rather 
strong. Understanding the origin of these constant density 
dark matter cores could be one of the most exciting chal- 
lenges facing astronomy in the next few years. 
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APPENDIX A: FULL PHASE SPACE 
ANALYSIS: THE TECHNIQUE 

Al The surface brightness profile 

For any population defined by a Michie-King model 
{Rh,nt,P) in a dark halo characteristic radius ro, we can 
define the quantity 



, ^ [ES^E™(i?,/ro;i?h,nt,^)-E°(i?0l 

xi = E - — 



Var[E°(J?, 



(Al) 



Here, the collection of points E°(J?i)) is the observed 
surface brightness profile, each of the points having its vari- 
ance Var[E°(i?,)]- The function E'"(f ; Eh, rit, /3) is the di- 
mensionless surface brightness of the Michie-King model 
with parameters {Rh,nt,P) at the dimensionless point f 
(measured in units of the halo scalelength) . This surface 
brightness is made comparable with the data by the dimen- 
sional coefficient E™, which is determined in such a way to 
give the smallest Xs- 



dxl/dE^ = 



(A2) 



We associate the quantity to the standard likelihood 



LE(Jih,nt,/3,ro) = exp (-xe/2) 



A2 The velocity dispersion profile 



(A3) 



Similarly, for any triple {Rh,nt,l3) and any pair of dimen- 
sional scales ro, po, we define the quantity 

2 _ [V^^^M^/ro■,Rl.,ntJ)~al,,%R^)f 

" Y VarKs°(a)] ■ ''^^^ 

Here, the collection of the points {Ri,aios°iRi)) is the ob- 
served projected velocity dispersion profile, each of the 
points having its variance Var[crios°(-Ri)]- The function 
S']2s{^j Rh,nt, is the dimensionless projected velocity dis- 
persion profile of the Michie-King model with parameters 
{Rh,nt,P) at the dimensionless point f (measured in units 
of the halo scalelength). This velocity dispersion is made 
comparable with the data by the square root of the poten- 
tial term $o, eqn ((7|, and in which the dimensional scales ro 
and po enter. We associate the quantity xi to the standard 
likelihood 

LCT(^h,nt,/3,ro,po) = exp (-x^/2) . (A5) 

A3 The tidal radius of the metal-poor stellar 
component 



As explained in Section [3.21 the two stellar components are 
allowed to have different tidal radii. Recall that we interpret 
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the tidal radius of the metal-rich as an approximate cut in 
energy of its phase space distribution. On the other hand, 
the tidal radius of the metal-poor stellar component is the 
true tidal radius of the dSph within the gravitational poten- 
tial of the Galaxy. At this radius, the density distribution 
of the halo is truncated too, according to the prescription of 
eqn. ([2]). Hence, such a tidal radius has to satisfy, at least 
approximately, the Roche criterion: 

Mdsph(r-t) _ Mmw{D - n) 



(A6) 



where MdSph('") is the mass profile of the dwarf, Mmw('") is 
the mass profile of the Galaxy and D is their separation. 

Suppose we have fixed a Michie-King model {Rh,nt,l3) 
and the pair of dimensional scales ro, po- Then the tidal 
radius of the model is 

rt'"(-Rh,Wt,/3,ro) = nt • f™(_Rh,/3) • r-o . (A7) 

We denote this by rip(_Rh, nt, /3, ro, po) and quantify the dis- 
crepancy using 

2 _ ('•r - ''t^) 

This compares the two different determinations of the tidal 
radius r™ and r^. The uncertainty Sr^ takes into account 
the uncertainty on the mass profile of the Galaxy, for which 
we allow the interval 2 > Mmw /{W^'^ Mq) > 0.5. We finally 
associate the quantity Xrt to the standard likelihood via 

Lrt(-Rh,nt,/3,ro,po) = exp (-Xrt/2) . (A9) 
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